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ABSTRACT 



Motivated by the observations of high Faraday rotation measures measured in 
cooling flow clusters we propose a model relevant to plasmas with comparable thermal 
and magnetic pressures. Magnetic field reconnection may play a major role in changing 
the topology of the magnetic field in the central cooling flow regions. The effect 
of the topology change is that cool flux loops can reconnect to hot flux loops that 
are connected to the overall thermal reservoir of the cluster. There can be a rapid 
recycling of mass between hot and cold phases on a time scale of ~ 3 x 10 8 - 10 9 yr 
which may reduce the inferred inflow and mass condensation rates by at least an 
order of magnitude. A central multiphase medium is a direct consequence of such a 
model. Throughout the cooling flow the filling factor of the hot loops (T > 2 x 10 7 K) 
is of order unity. The filling factor of the cool loops(< 2 x 10 7 K) is 0.1 — 1% with 
a corresponding mass fraction of cold phase of 1 — 10%. A crucial parameter is the 
coherence length of the field relative to the cooling radius and the distribution of field 
energy with scale. When the cooling radius is greater than the field coherence length 
then cooling flows proceed as usual. When the coherence length is greater than the 
central cooling radius, however, the thermal energy of the reservoir can be tapped 
and the mass condensation rates may be very significantly reduced. Three additional 
conditions must be satisfied: I. Cold loops must be able to fall at least as far as the 
mean distance between hot loops in a cooling time; II. Loops must enter an evaporative 
phase on reconnecting; and III. A sufficient number of hot loops penetrate the cold 
phase region to power the radiative losses. 
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1. INTRODUCTION 



Calculation of the mass flow rate into the centers of clusters of galaxies with short cooling 
times gives large values of order ~ 100 — 1000M Q yr _1 (Binney and Cowie 1981). This point has 
dominated most of the thinking on the thermal evolution and central gas dynamics of clusters of 
galaxies (Fabian 1994). We propose a mechanism that will reduce these mass inflow rates by an 
order of magnitude in clusters where the coherence length of the field is greater than the cooling 
radius. The mechanism we invoke here involves, however, the study of turbulent magnetized fluids. 
Given our state of knowledge of the properties of these flows our calculations can, at present, only 
be regarded as indicative. However, the stabilization does seem to be a reasonable mechanism and 
we present it here in the spirit of the pioneering feedback work of Tabor and Binney (1993) and 
Binney and Tabor (1995). 

Our motivation for rethinking this subject was the beautiful data from Ge (1991), Ge and 
Owens (1993), Owens et al. (1990) clearly showing that cooling flow clusters have rotation 
measures between 10-100 times those that have no established cooling. The inferred field 
strengths are at or near equipartition values with the thermal pressure. Magnetized clusters 
have been studied before in detail by Soker and Sarazin (1990). They assumed that once the 
magnetic field built up it was annihilated down to equipartition values via reconnection process. 
Normally, however, reconnection does not annihilate magnetic fields (Taylor 1974). Even with fast 
reconnection, annihilation of significant volumes of magnetic energy density does not take place. 
To dissipate the field globally requires a global anomalous resistivity and this is very difficult both 
to achieve and sustain except in very small volumes. 

The issues we have to address are complex. Magnetic field reconnection normally rearranges 
the field and may allow some change in the mass-to-flux ratio of loops. Since the magnitude of the 
field does not change significantly during reconnection the magnetic pressure is not significantly 
altered. The crucial point of physics for this work, however, is that reconnection significantly 
alters the topology and connectedness of the magnetic field. 
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In the central regions of cooling flows the line widths are random (Heckman et al. 1989) and 
it is most likely that the fluid is turbulent. The most recent analysis of the nature of thermal 
conduction in cooling flows (Tao 1996) indicates that there may be some but not complete 
suppression of the thermal conductivity. Additional suppression may result from electromagnetic 
instabilities in the cluster gas (Pistinner et al. 1996). We will adopt the values consistent with 
Tao's work in our subsequent calculations although we regard the thermal conductivity in a 
tangled field as still an unsettled question (cf. Rosner and Tucker 1989, Tribble 1989, Chun and 
Rosner 1993). 

Recent spectral X-ray measurements have established that the X-ray emitting gas is composed 
of two (or more) phases (Fukazawa et al. 1994; Fabian et al. 1994). We incorporate this important 
new datum into our model. Our proposed solution to the suppression of the large mass drop out 
rates in cooling flows is to recycle the cooling ( < 10 7 K) gas back into the hot phase (> 10 7 K) 
and then allow it to cool again. Only a small fraction of the total cooling mass need be in the 
cool phase to match the observed low-temperature X-ray emission. The mass condensation rate 
is, then, the amount of cool gas per recycling time. Since current estimates of the amount of 
condensed mass are calculated assuming that the cooling gas eventually leaves the flow with no 
recycling, the currently estimated masses can be suppressed by the ratio of the recycling time to 
the age of the cluster. 

More specifically, the recycling solution we propose is as follows: 

The magnetic field lines engulfed by the outward moving cooling wave have a spectrum of 
scales. The ones with scales greater than the scale of the cooling wave remain connected to the 
essentially infinite heat reservoir of the cluster. Even with the turbulent thermal conductivity 
given by Tao (1996), field lines topologically connected to the outside of the cooling wave region 
remain hot and can be treated as isothermal. These hot tubes are the pipes that the cooling flux 
loops can tap into. Smaller scale cooling loops are not connected to the heat bath of the cluster. 
However, small-scale cooling loops can be heated and recycled by the process of magnetic field 
reconnection giving the required recycling mechanism. 
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In the next section (2) we set up a simple physical model and present estimates of the 
characteristics of our solution. Thermal structure and recycling are discussed in more detail in 
section 3. In Section 4 we discuss some observational consequences and a summary and conclusions 
are given. 



2. BASIC MODEL AND ESTIMATES 

Assume an initially tangled isotropic field across the cluster with a plasma parameter 
P = (-fthcrmal/-fmagnetic ) = A) > 10 since the cluster is initially dominated by thermal pressure 
and not magnetic pressure. We assume there is a spectrum of loops as a function of scale. We 
envisage these intracluster fields to be generated by radio galaxies, outflows from AGNs and 
starburst galaxies and the outflow processes from normal galaxies that enrich the intracluster 
medium that are presumably from a supernovae driven source. With these processes in mind we 
estimate that the scales of the fields in the general intracluster medium can range from ~ 10 kpc 
up to ~ 1 Mpc. The field may be completely coherent across the cluster if there is a large scale 
primordial field or if large scale fields are generated in the process of structure formation (Kulsrud 
1995). 

After a time t the material enveloped by the outwardly propagating cooling wave (with radius 
defined by that radius at which i coo i = t) can cool and move inward (Bertschinger 1989). In the 
kinematic inflow phase the fields are frozen in and the stretching and amplification that occur are 
well described in the paper by Soker and Sarazin (1990). The field becomes increasingly radial 
and we will assume for simplicity of analysis that it is exactly radial in this phase. 

Assuming complete ionization, the cooling law that describes the propagation of the cooling 
wave is given by: 

3 tin hT 

n A(T), (1) 



2 dt 

where T is the gas temperature, n is the density, and A(T) is the cooling function. We assume 
that the gas and dark matter both have the distribution of an isothermal sphere with a gas 
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temperature, T g , and dark matter temperature, Tj, with gas density 

n = n ( r f) a . (2) 

where a = 2T c i/T g and the subscript refers to gas properties at a fiducial radius r . The 
outwardly propagating cooling wave has the solution: 

( * 

r coo \ = r , (3) 

\lcool,0 / 

with 

/ r ^ ' 



^cool = Ucool,0 ( ) , (4) 

V r cool / 



where 

3kT 



^cool,0 



(5) 



_2n A(T 

and f coo i,o = ( r o/*cooi,o)- The mass flux into the cooling bubble, at any given time t and at radius 

^coob ls: 

Mcooi = M (^) 3 " 2a , (6) 

where 

M =(^]. (7) 

\ J cool,0 / 

As the cooling wave propagates outward, the material trapped in small magnetic loops will start 
to cool and fall inward. 

The plasma parameter (3 varies as a function of radius as: 

a/3 

TO 

and the flow can be significantly magnetized when (3 ~ 1. We are assuming that the flow is 
sufficiently large as to entrain field lines into the central regions of the cooling flow, allowing a 
build up of the magnetic pressure. The increased magnetic pressure does not change the estimates 
of the accretion rates that are based purely on the thermal properties, not the momentum balance. 

As discussed by Bertschinger and Meiksin (1986) and Meiksin (1988), if the thermal 
conductivity is within an order of magnitude of the Spitzer value, the cooling along flux tubes that 



P = Po ( - ) , (8) 
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are simply connected to the outside of the cooling wave is negligible. Therefore these connected 
flux tubes do not contain any cool material. The flux tubes that are not connected can cool. 
This leads to a two-phase medium of hot connected flux tubes and cool disconnected tubes. As 
discussed by Soker and Sarazin (1990) even if all the flux were reconnected and dissipated the 
associated heat input would have negligible effect on the thermal balance. 

Reconnection is a natural way to recycle the cold material back into the hot phase. The exact 
details are unclear but the result is probably bounded by the rules for collisions of flux tubes (cf. 
Norman et al. 1996) as a lower bound or explosive reconnection with all the mass injected into the 
background medium between hot and cold phases as an upper bound. 

Let us assume that a fraction, of the mass and flux is in hot tubes and the rest in cold 
tubes. Let us also assume that there is a mean flux reconnection rate, T, giving a mean mass 
exchange rate M e between hot and cold of: 

M e = (^] 'r*, (10) 



M. 

where the flux to mass ratio must be taken in mean averaged sense here. Therefore since the 
timescale for mass exchange is T -1 , the fraction of the total mass, Mh, within a cooling radius 
that is in the form of cool condenstates ( < 10 7 K) is: 

M rec = (Try 1 M h , (11) 

where r is the time for mass to leave the hot phase, which is on the order of the age of the cluster 
because of thermal heat conduction. The amount of mass being recycled at any given time is 
negligible compared to the total amount of mass in the hot phase. 

Concentrating on one loop and using the fast reconnection law, Soker and Sarazin (1990) give 
a typical reconnection timescale, tr = T~ l ~ l t /{eVA) of: 



r*=iov^-j 



s i u w u . / „ . (12) 



cm 3 



i 



for our parameters. We assume that reconnection occurs with a reconnection velocity of e times 
the Alfven velocity. The characteristic scale length for reconnection is here taken to be the 
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transverse dimesion of a loop. Note that when large scale hot loops reconnect with small-scale 
cool loops the reconnection time given in (12) should be regarded as an upper limit since the 
reconnection time may be considerably shorter due to the smaller size of the loops. The internal 
Alfven speed Va in the cool and hot loops does not change significantly since Va oc n 7 where 
< 7 < 1/6. Complete reconnection will occur on a single collision if the reconnection time is 
less than or equal to the crossing time of the loops. If this is not satisfied then collisions with hot 
loops will continue until the process is complete. This is similar to shredding. For more details 
regarding incomplete reconnection during flux loop collisions see Norman et al. (1996). Generally, 
the complex process of shredding of loops and blobs as they fall into the central regions of the 
cooling flow is one where cool gas distribution cascades to smaller and smaller scales. A small 
inefficiency in the recycling would permit some material to continue to cool and make stars, but 
at a much reduced rate compared to the standard cooling flow prediction. 

We note here that buoyancy is an inefficient means for transporting flux out of the cooling 
region (Soker and Sarazin 1990). Flux loops will be buoyant and reach an appropriate scale 
height. We assume that most of the loops are still trapped, although there is probably still some 
flux loss. Reconnection will then occur and the value of mass and flux for the product loops will 
change. There will then be a rising and falling of loops as the reconnection process continues 
and re-sorts the mass and flux, although it is not possible to follow an individual piece of flux 
and mass. Magnetic tension in reconnected loops may lift the cool optically-emitting material 
out of the center of a cooling flow, accounting for the wide range of scales over which the optical 
filaments are observed (Zoabi et al. 1996). In practice, loops will shred into smaller structures due 
to Kelvin-Helmholtz instabilities after only a few Brunt- Vasaila oscillation periods (Reale et al 
1991). As discussed below this will enhance their evaporation and mass recycling. 



3. THERMAL STRUCTURE AND RECYCLING 
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We now describe the thermal structure of the intracluster gas. The radiative cooling of the 
hot phase is nearly balanced by heat conduction from the outer portions of the cluster into the 
cooling region by those extended flux tubes that connect the external heat reservoir to the cooling 
interior. While perfect balance between cooling and heat conduction requires too high an external 
temperature compared to observations, acceptable solutions are found when a small amount of 
inflow is allowed (Bertschinger & Meiksin 1986), provided the temperature decline is not more 
than by a factor of ~ 2 — 3 (Meiksin 1988). Accordingly, we may expand the fluid variables in the 
small parameter e = Mfl ow /M^, where Mfl ow is t ne t rue accretion rate, and M^, given by equation 
(7) above and evaluated at the cooling radius, is the rate that would be inferred from the X-ray 
luminosity assuming heat conduction were suppressed. Thus, for example, we may write, to first 
order, n(r) = no(r /r coo i) _a [l + eni(r)] and T^ir) = To(r)[l + eT\(r)]. The zeroth order equation 
for pressure support (both thermal and magnetic) then requires, for (3 > 1, 



OL = T^, — 



IT v i —1 

-'"l+Pir)- 1 ] , (13a) 



when the gas temperature varies slowly with radius compared to the density. We point out that 
in the case ~ and where the cooling flow is sufficiently magnetized with [3 ~ 1, then a ~ 1. 
Consequently, from equation (6), M C ooi oc r, a case favored by the observations. We will regard a 
below as being directly determined from the measurements of the X-ray surface brightness of the 
hot component, since the radial dependence of the surface brightness depends only weakly on that 
of the temperature. 

Setting x = r/r coo \ and y = To(r)/Th, the zeroth order thermal energy equation becomes 

where Ap, c = ^(T/jT/i/noA)] 1 / 2 is the Field length, evaluated at the cooling radius (Begelman 
and McKee 1990; Field 1965). Near T ~ 2 x 10 K, the cooling function is nearly independent of 
T. The Field length is then 

0.003 cm" 3 \ / T h \ 7 / 4 



n 



A F , c ^100kpc(— ) (aTlJjTK) ' ( 14 ) 
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where we have used the cooling function of McKee and Cowie (1977). The solution of equation 
(13b) is then 



7 / '"cool 



T(r) = T h 
for a ^ 3/2 and a^l; 

T(r) = T h 

for a = 3/2, and 



1 



4 V A F , C J (2a-3)(a-l) W a ~ 2 x) x y °° \ x 



2/7 



(15a) 



7 ( r CO oi \ logz 1 
2 I Af,c / x x 



+ ?/co /2 (l 



2/7 



(156) 



T(r) = T A 

for a = 1. Here, is an integration constant 



2/7 



(15c) 



The solution given in equation (15) is a gradually decreasing temperature profile toward the 
center. A typical cooling flow will have a cooling radius of ~ 100 kpc and Tj,~5x 10 7 K, giving 
^cooi/Af, c ~ 0.3. For the simple, observationally favored case of a w 1 and ps 1, we find 
that the temperature at r ~ 0.1r coo i will be below its value at the cooling radius by a factor of 
~ (1 — 0.3 log 10) 2 / 7 « 0.7, consistent with ^(SCA measurements for the hot phase (Fukazawa et 
al. 1994). 

A strong argument in favor of high mass drop-out rates is the high luminosities measured in 
the low-temperature X-ray emission lines, such as the Fe L lines (Mushotzky et al. 1981; Fabian et 
al. 1994). These lines demonstrate the presence of low temperature gas that is cooling at a high 
rate, and buttress the argument for a correspondingly high rate of mass condensation out of the 
flow. We next show how recycling may significantly alter this point of view. 

Consider two phases in pressure equilibrium, a hot phase and a cold phase at respective 
temperatures Th and T c , with number densities nh and n c . The cold phase arises from thermally 
unstable gas suspended in magnetic flux loops. Let f c denote the volume filling factor of the cold 
phase. The mass M c in cold phase is related to that in the hot, M^, by 
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M c = -±-^M h *f c ^M h . (16) 
1 - f c n h T c 

The mass cooling rate M c of the gas in the cold phase is then related to that in the hot, Mh, by 

for a cooling function A(T). In principle, it is possible for the cooling rate inferred from emission 
lines generated by the cold phase to exceed that inferred from the hot. We will show below, 
however, that in our model these rates should be comparable, though some variance in the ratio is 
to be expected. 

The crucial point now is that as the flux loops carrying the cold phase sink they will reconnect 
to the hot flux tubes. Thermal heat conduction carried along the hot flux tubes will then quickly 
reheat the cold phase gas to Th, resulting in a continual recycling of the cold phase material. The 
actual mass drop-out rate will thus be much less than inferred from the line luminosities. This is 
possible provided three criteria are satisfied: 

I. The cold loops are able to fall at least as far as the mean distance between hot flux 
tubes in a cooling time: i>//t C ool,c > ^h] 

II: Loops must enter an evaporative phase on reconnecting; 

III. There is a sufficient number of hot flux tubes penetrating the cool phase region to conduct 
heat inward from the infinite heat bath surrounding the cool phase to power the total radiative 
losses from both the cool and hot phases. 

These conditions may be met as follows: 

I. To make a simple estimate assume that for the purposes of argument that the cooling 
time is of the order of the time for a loop to cross the cooling region. Then for a volume filling 
factor of the hot phase within the cooling region of fh ~ 1 , the covering factor of the hot phase is 

2 

~ fh ~ 1- The cool loops will then collide with the hot phase on the order of a cooling time of 
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the cool phase. 

II. The cool loops must have a characteristic scale less than the Field length when they 
reconnect to the hot loops. For loops that are initially magnetically isolated and collide early 
on ( when they have temperatures of order, say, 3 x 10 6 K) with hot loops they will undergo 
evaporation in the classical limit and be recycled since 

. /0.03cm- 3 \/ T h \ 5 / 2 / T c \ L6 

AF = 3kpC {-^^)[3^WK) [J^Wk) ' (18) 

where we have allowed for the radiation of the cooler material using the cooling approximation 
given in McKee and Cowie (1977). 

Evidence for absorption by very cold clouds has been argued for on the basis of the X-ray 
spectra in some cooling flows, with inferred characteristic cloud sizes of ^ 1 pc (Fabian 1994). In 
our model, such clouds will ultimately result from any small inefficiency in reheating. We suggest 
that these clouds are shortlived, and will evaporate when they reconnect with hotter loops. For 
the small cloud sizes suggested, heat conduction occurs in the saturated limit, so that the clouds 
will disappear in an evaporative wind (McKee and Cowie 1977). The timescale for the clouds to 
evaporate is 

where R is the cloud radius, and the saturation parameter, ao = 4(7^/3 x 10 7 K) 2 /(n/ l i? pc s ), is 
~ 10 3 , taking 4> s ~ 0.3 (Balbus and McKee 1982). 

From analyses of absorption-line studies in both X-ray and radio bands the temperature of 
the clouds must be less than 10 K (Crawford and Fabian 1992; Dwarakanath et al. 1994). The 
X-ray absorption measurements indicate H I column densities of the order of 10 20 — 10 21 cm -2 
(Fabian 1994). At the canonical pressure of the central cooling flow this corresponds to scales of 
only 10 16 — 10 17 cm corresponding to cloud evaporation timescales of ~ 10 5 — 10 6 yr. Therefore, 
when the reconnection occurs the cold H I will evaporate and be recycled into the hot phase. 
The cloud destruction time is given by the time for reconnection with hot loops (~ 10 8 yr). For 
unit covering factor the cloud destruction rate in the central region of order -R C entrai ~ 1 kpc is 
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~ l(i?ccntrai/lkpc) 2 (iVHi/10 21 cm- 2 )M yr" 1 . 

III. The hot phase will occupy most of the volume. A reasonable estimate of the cooling 
filling factor is f c « (rV) _1 (T c /T h ) ~ 1(T 3 << 1. From equation (16), only about 1% of the total 
intracluster gas within the cooling radius will then be in the cool component. Because of the 
relatively greater efficiency for generating X-ray photons per unit mass of the cool phase compared 
to the hot, however, equation (17) shows that the mass cooling rate of the cool component will be 
comparable to that of the hot, and may even somewhat exceed it. If the cool phase were detected 
through an X-ray emission line, such a disagreement, which could be interpreted in the standard 
cooling flow picture as an abundance anomaly of the element producing the line, could equally 
well be interpreted as evidence for a second, denser and cooler phase. 

4. SUMMARY, CONCLUSIONS, AND OBSERVATIONAL CONSEQUENCES 

We have proposed a model for significantly reducing the mass inflow rates in clusters of 
galaxies. Reconnection processes acting on the topology of the field regulates the inflow of heat 
from the vast heat bath of the entire cluster. Reconnection has no energetic advantage at all; it 
is simply the topology change that allows heating of the matter in cold loops and a relatively 
fast recycling time. Therefore the cool gas can be used many times, in contrast to conventional 
cooling flow models in which once the gas cools, it necessarily becomes invisible and can never be 
recycled. A crucial parameter is the coherence length of the field relative to the cooling radius 
and the distribution of field energy with scale. When the cooling radius is greater than the field 
coherence length, cooling flows proceed as usual. When the coherence length is greater than the 
central cooling radius, however, the thermal energy of the reservoir can be tapped and the mass 
drop out rates can be very significantly reduced. Three additional conditions must be satisfied: 
I. Cold loops must be able to fall at least as far as the mean distance between hot loops in a 
cooling time; II. Loops must enter an evaporative phase on reconnecting; and III. A sufficient 
number of hot loops must penetrate the cold phase region to power the radiative losses. 
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Central regions will be multi-phase with a hot-phase filling factor of order unity and the 
cold-phase filling factor very much less than unity. There is rapid mass circulation on a timescale 
of order 10 8 yr. 

Filaments will be highly magnetized. In general, the central regions of cooling flows should 
look very magnetized and filamentary. Detailed rotation measure maps could map the central field 
geometry. The two-phase model we outline here is consistent with the recent ASCA observations 
of the centers of cooling flows (Fabian et al. 1994) although many details remain to be understood. 

We thank our colleagues for many useful and enjoyable discussions on this subject in 
particular F. Adams. S. Cowley, M. Donahue, T. Heckman, D. Melrose, R. Mushotsky, R. Sudan 
and M. Voit. We thank the Institute for Theoretical Physics at UCSB for the stimulating 
environment where this interdisciplinary effort was possible. The ITP is funded by NSF Grant No. 
PHY94-07194. A.M. is grateful to the W. Gaertner fund at the University of Chicago for support. 
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